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1 Introduction 



In a previous paper |CaKi2010] , the author and Holly Carley developed a constructive 
approach to the prescribed mean-curvature equation |GiTrl983] 



for hypersurfaces which are graphs of some scalar function u over M", having mean- 
curvature H G (Cg nL^)(M^) suitably small. In ([T]) the + sign refers to hypersurfaces 
in Euclidean M"+^, while the — sign refers to spacelike hypersurfaces in Minkowski 
spacetime M^'". This latter setting, when n = 3, is equivalent to the electrostatic 
Maxwell-Born-Infeld problem for the electric potential (p (x u with regular charge 
density p (x H and small Born constant /3. In the present paper this constructive 
approach is generalized to the stationary Maxwell-Born-Infeld field equations for 
the electric and magnetic potentials and A, or rather their electric and magnetic 
fields E and B, with regular charge density p and current vector-density j. We here 
are interested only in electromagnetic fields in unbounded space with vanishing 
conditions at spatial infinity. 

Since the Maxwell-Born-Infeld field equations were proposed as a nonlinear rem- 
edy for the infinite field energies and momenta incurred with point charge sources in 
the linear Maxwell-Maxwell field equations (the Maxwell-Lorentz theory), our study 
of the Maxwell-Born-Infeld field equations with regular source terms p and j may 
seem irrelevant from a fundamental physical perspective. However, sufficiently far 
away from delta function sources the fields should be practically indistinguishable 
from regularizations thereof, yet very little is known about solving the Maxwell- 
Born-Infeld field equations for non-special choices of source terms. 

In the next section we briefiy summarize the stationary Maxwell-Maxwell and 
Maxwell-Born-Infeld field equations in dimensionless units in which the speed of 
light c = 1. There we also explain our terminology. Section 3 puts forward our 
perturbative strategy to solve these equations. In section 4 we prove the convergence 
of the perturbative, formal power series solution. In section 5 we specialize to the 
electrostatic and magnetostatic subproblems. Section 6 concludes the main part of 
this paper with a brief summary and an outlook on applications. The appendix 
summarizes some important details concerning the recurrence relation for the coef- 
ficients of the perturbative power series which were supplied by my colleague Doron 
Zeilberger in collaboration with Shalosh B. Ekhad HI. 
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2 Stationary electromagnetic Maxwell— Born— Infeld theory 

Maxwell's classical electromagnetic theory comprises Coulomb's law, Faraday's law, 
the Ampere-Maxwell law, and the unnamed law of the absence of magnetic charges. 
In stationary situations the Ampere-Maxwell law reduces to Ampere's law, and Fara- 
day's law to its stationary special case. Coulomb's law states that an electrostatic 
charge density p in M? is the source for the electric displacement field D, 

V ■ D = ATxp Coulomb's law , (2) 
while Faraday's law says that the electric field when stationary, is curl-free: 

V X = Faraday's law (stationary) . (3) 

Similarly, Ampere's law states that a solenoidal electric current vector-density j in 
is the source for the magnetic excitation field H, 

V X H = Airj Ampere's law , (4) 
while the law of the absence of sources of the magnetic induction field B reads 

V ■ -B = no — magnetic — charges law . (5) 

The four fields E, D; B, H need to be linked by an "aether law," with the help of 
which one can eliminate one of the electric-, and one of the magnetic- type fields from 
the equations and obtain a closed, nonlinear set of first-order pde for the two remain- 
ing fields. Note that for us "aether" is just convenient shorthand for "electromagnetic 
vacuum;" of course, for Maxwell and his contemporaries it had a substantial meaning 
until it was demolished by Einstein. Maxwell's law of the "pure aether" reads 

E = D Maxwell's E - law; (6) 

B = H Maxwell's B - law; (7) 
whereas the aether law of Born and Infeld |BoInl933/34| says 
jj B'^D ■ HH 

E = , Born - Infeld's E - law; (8) 

v/1 + fiWE)\^ - - ■ Df 

B = ^ Born - Infeld's B - law. (9) 

v/1 + PWD\^ - \H\^) - f3^{H ■ Df 
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In ([H]) and (3 G (0, oo) is a hypothetical new constant of nature (in the dimen- 
sionless notation of |Kie2004a] ). which we have called Born's aether constant. In the 
limit /3 — )■ the Born-Infeld law goes over into Maxwell's aether law. 

To avoid any ambiguity as to which set of field equations we refer to in this paper 
(and elsewhere) we will always use the combination "Maxwell- whose aether law." 
Thus, in particular, we will speak of the Maxwell-Maxwell equations when normally 
one just refers to the Maxwell equations. 



Assume that the sources p G (Cq fl L-^)(M^), and j G (Cg H L-'^)(M^) componentwise. 
We seek stationary solutions in unbounded space of the electromagnetic Maxwell- 
Born-Infeld field equations satisfying vanishing conditions at spatial infinity. 

We begin by recalling the just mentioned fact that for /3 — )■ the system of 
Maxwell-Born-Infeld field equations reduces to the system of Maxwell-Maxwell field 
equations, which for the stipulated class of sources and asymptotic conditions are 
uniquely solved by Ec = Dq and = Ha, where 



Here, the subscripts c and a stand for Coulomb and Ampere respectively. This 
suggests that for /3 > we may seek a solution to the Maxwell-Born-Infeld field 
equations for the same sources and asymptotic conditions as for the Maxwell-Maxwell 
field equations by setting Dbi = Dc + Dp and Hbi = Ha + Hjs, with lim/3_j.o -D/3 = 
and lim/3_s.o = 0, from which the electric field Ebi and magnetic induction field 
Bbi are obtained through the Born-Infeld aether law. 

To determine and if^ we first note that, since V • Dc = 47rp, we have that 



3 Formal perturbation theory 




(10) 




(11) 



V-Dis = 0, 
and since V x Ha = 4:7^ j, we have that 



(12) 



VxH, 



P - 



0. 



(13) 



Thus, is a curl field, while Hp is a gradient field. 
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Incidentally, since moreover Dq is a gradient field and H\ a curl field, we now 
see that the Ansatz Dbi = Dq + Dp and H^i = Ha + Hp corresponds precisely to 
the Helmholtz decomposition for the fields D^i and ifei- For each field we have one 
of the two Helmholtz components and seek the remaining one. 

Coming back to the determination of Dp and Hp, we now need a curl equation 
for Dp and a divergence equation for Hp. So, using (|H]) to eliminate E in favor of D 
and if in ([3]), and using ([9]) to eliminate B in favor of D and if in ([5]), we obtain 
a system of nonlinear vector pde for Dp and Hp which is closed conditioned on Dc 
and Ha being given. This systems of conditional vector pde reads 

Vx ^ D'^*D-HH ^ 

+ /3^(|i^P - \H\^) - P\H ■ DY 
H + (5^H -DD 

V- , ^ =0, 15 

v/1 + ^^{W - W) - P^{H ■ Df 

where D and H stand for Dc + Dp and H for Ha + Hp, respectively. 

The analyticity properties of the nonlinearities now suggest in particular that a 
solution pair Dp, Hp for the nonlinear first-order vector problem for small but finite 
(3 is itself analytic in /3, i.e. satisfies a power series Ansatz, more specifically 

oo 

Dp{x) = Y,P"'D''''\^)^ (16) 
fc=i 

oo 

Hp{x) = J2^"'H('\^)^ (17) 

k=l 

with each D^''^ and if'-^^ independent of (3. Inserting this power series Ansatz into the 
pair of equations (IT^ . ( ITS!) and sorting according to powers of (3, we find a hierarchy 
of linear equations, the k-th members of which read: 

k 

V X D'^^^ = V X ^ (^X^h) j^(k-h) ^ y (/^)^(fc-/^) j . ^ ^ |^ (^g) 
h=l 

(together with V ■ D^'''> = 0), and 

k 
h=l 

(together with V x H^^^ = 0); here, D^°^ := Dq and ii(°) := Ha, and X^''^ and Y^^^ 
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are the coefficients of ^^'^ in the power series expansions of the two formal functions 

^ -1/ + - 1^1') - ■ Df 

and 



respectively, with D = Dq + Dp and = + Explicitly, and including /i = 0, 
we have = —1, while for /i G N we have 

= - ^ M,^(-l)^-' A') (20) 

with 

2i 

|6|4j=/i-2in=l 

i 

and, for < Z < j, 

I 20-0 
^ JJ(L»-//)("2m-l).(^^^)(a2^) JJ 2?(62n-l).^(62n)_ (23) 

I''l2i + l''l4(i-i) "^=1 "=1 
=h+i-2j 

Here, \1\k '■= X^^i, with taking any non-negative integer values; so in particular, 

i=l 

^\e\K<o - = 0. Also, Mo := 1 and 
(27 - 1)!! 

Mj^{-iy^-^^^; jeN (24) 
are the Maclaurin coefficients for + z. Having X'-^^ for /i = 0, 1, 2, we find 

= - J^X^^-^-^) J] L>('^i^-//("2); /leN. (25) 

5=0 |a|2=ff 

Note that X^'^^ and y^'*) contain only terms of order < h oi |-D^P, and Dp- Hp. 
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The expressions for X^'^^ and F*^*^^ look somewhat unwieldy, but since there is 
scant empirical evidence (if any!) for the need of corrections to the solutions of the 
Maxwell-Maxwell field equations in the regime of classical phenomena, for many 
practical applications it may suffice to just compute the first correction term to 
the Coulomb and Ampere fields. These are of a more manageable structure. In 
particular, D*^^^ satisfies 

V X D« = V X (|(|/^c|' - \Hk\^)Dc + Dc ■ H^Hj,^ , (26) 
together with V ■ D^^^ = 0, and H^^^ satisfies 

V ■ if« = V ■ (i(|/^c|' - \Ha\^)Ha - Dc ■ HaDc) , (27) 

together with V x if^^^ = 0. 

As to the solvability of the linear first-order PDE flTSj) . f lT9|) . assuming that their 
right-hand sides are in Cg , this pair of PDE has a unique solution in Cq'" given by 

k 

D(k) = pJ2 (^x^^'^d'^^-^^ + j , ken (28) 

h=l 
k 

H^^) = Q ^ (^x^h)jj(k~h) _ jjik-h)^^ ^ ^ ^ (29) 
h=l 

where P : Cq'" — )■ Cq'° projects onto the solenoidal subspace of Cq'", and where 
Q : Cq'" — Cg'" projects onto the curl-free subspace of Cq'". More explicitly, for a 
7o''" with V ■ V^^^ e C°'' 



vector field G Cn'" with V ■ V^^'^ G Cn'" n L\ we have 



PV^^\x) = V^''\x) + V /" ' d^y, (30) 

J \x-y\ 

and when V x V'^''^ G (7°'" n we have 

QV^'\x) = V^''\x) - V X / ^^ ""^^'^^y^ d^y. (31) 

J \x-y\ 

To summarize, so far we have seen that the Ansatz D = Dq + X^fc^i and 
H = Ha + Yl=iP^^H^^\ with D^^^ and H^^'^ for A; G N recursively given by ([28]) 
and f l29|) . yields a formal series solution pair for the stationary Maxwell-Born-Infeld 
field equations with regular sources. We next prove its convergence in Cq'" for small 
enough /3, given the sources p G Cq and j E Cq. 
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4 Rigorous perturbation theory 



modulo some MAPLE calculations 



4.1 Cq'° consistency of the formal perturbation series solution 

We first confirm that the formal series D^s = Yl'k=i l^'^^ ^^'^^ ^^'^ ~ Yl'k=i /^'^'^^^''^ 
are Cq'"^ series. For this we have to check that all D*^^^ and H^''^ as given in flHU]) and 
( 13T|) are in Cg'". But this follows inductively from the facts that D^^^ = Dq G Cq'" 
and = Ha G Cq'", that Cq'" is a Banach algebra, and that P : Cq'° — Cq'" and 
Q : Cg'" — )■ Cg'"; hence, each partial sum of Djs = /S'^^D'^''^ is in Cg'" and so is 

each partial sum of if/3 = Yl'kLi P^'^H'^''\ 

4.2 Absolute convergence of the formal perturbation series solution 

Since ||D^|| < E^i /^^'^P^'^^ll and \\H,3\\ < E^i /S^ll^^'^^ll in the sense of partial 
sums, to prove absolute convergence of the power series for Dp and Hp it suffices to 
show that + ll-^*''^''!!) < sufficiently small /3, given p and j. 

We first estimate \\D^'''^\\ + \\H^'''^\\ for A; e N in terms of the \\D^^^\\ + for 
all I < k. Beginning with for G N we find 

^ IIP^''' ^x*^'^'/}*'*''"'''' + y'-'^-'iJ*'^"''-' j II (32) 

< ||5^'_^(xWZ}(^-'') + II (33) 

< ^'_J|(^x(^)z}{'=-'') + rWij{^-''))|| (34) 

- + "■^^''^^^''"^^0 ^^^^ 



Here, inequality f p3|) holds because P : Cg'" — >■ Cg'" is a projector, inequalities ([3i 
and 035]) are just the triangle inequality, while fl36p is valid in Banach algebras (here 
Cg'"). Similarly, for ||if'''^'*|| and all /c G N we find 

Adding these estimates for H-D^'^-'H and ||-?i^'^''|| yields 

\\n(k)\\ + ||/jW|| < 5^^_^(||xW|| + ||r(")||)(p('=-'^)|| + ||ii(^-'^)||). (38) 



Since ||-'^^'^^|| and depend in turn on ||-D*-'''|| and ||-f^'''''||, we need to carry on 

and estimate \\X(^^\\ + ||yW|| in terms of the + \\H^^ for all I < k. As to 

we have = 1 and 

ll^^'^ll<El^.lE(j)ll^/..vll (39) 
i=i 1=0 ^ ^ 

and, for h 

\\Zk,,o\\ < E nil^^'^"-^^llll^^''^"^ll (40) 

\b\4j=h-2j n=l 
4i 

|6|4j=/i-2in=l 

\\Zk,jA< E flll(^-^)^"""^^llll(^ + ^)^"'"^ll (42) 

|a|2j=fe-j"^=l 
2i 

^ E n(ii^^"'"^ii+ii^^"'"^ii)' (43) 

\a\2j=h-jm=l 

and if J > 1, which can happen only when h >2, then for < / < j we have 

\\Z^.,J,l\\ < E nil(^-^)^"'""^llll(^ + ^)^"''"^lin 11^^'"'-^^ II 11^^'^"^ II (44) 

=h+l-2j 

- E n (ii^^"""^!! + ii^^"'"^!!) n + ^^^^ 

=h+l-2j 

Finally as to ||F('*)||, we have 

h-i 

S=0 |a|2=g 
h-1 

< Ell^^''"'"''^ll E (ll^^^'^ll + ll^^"'^ll)(ll^^"'^ll + ll^^"'^ll)- (47) 

5=0 |a|2=fl 
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Having completed our estimate of H-D'-'^^ll + H-f^*-^^!! for A; G N in terms of the 
||/)(0|| + ||ij(0|| for all I <k,we now abbreviate \\D^'^'^\\ + = A^^") and recast 

our estimates more user-friendly as 



h-l 



-h) 



(48) 



h=l 9=0 

together with =1, then 

2 

|a|2=0 m=l 



|a|2=9 



(49) 



and, for h >2, 

\\X^'"^\\ < \Mi\ 



i=2 



|b|4=/i-2n=l |a|2=/i-l '"1=1 

|6|4j=h— 2j n=l |a|2j=/i— j m=l 



(50) 



film) 



i-1 



2/ 



40-0 



+E ^ E n'v'-'n 

/=1 ^ ^ |a|2i + IM4(j-i)"l=l "=1 
= h + ;-2j 

By an inductive argument we now estimate each A^*^'^^ in terms of the 2k + 1-th power 
of Ar(o) =: Af. 

In this vein, setting A; = 1 in (HH!) we obtain the estimate 

iV« < (llX^II +A/'')a/', (51) 

and for we have fj^9l) . so 

iV(i) < lAf\ (52) 
Next, suppose that for all k = 1, ...,k^, there exists some Rk such that 

N^^^ < RkAf^''+\ (53) 
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Also set Rq := 1 (for, N^^^ = J\f). Inserting these estimates into (HSj) then yields 

fc,+l h-l 

/i=l 9=0 \a\2=9 

As to the IIX^'^-'II, for h = I we already have the estimate (H9|) . and inserting (l53l 
into ( [50|) we find for /i > 2 that 

"4 2 " 



IIX^II < |Mi| 



J=2 



|fe|4=/i.-2n=l |a|2=h-lm=l 

|ti|4_,=/i— 2j ra=l \a\2j=h~j m=l 



(55) 



2am+l 



J-1 



2/ 4(j-/) 
\'^\2l + \''\4{j-l)^=^ "=1 



\a\2=h—l m=l 



\b\4=h-2 n=l 



(56) 



+Ei^^-i E n^^"+ E 11^- 

j=2 l\b\4j=h-2j n=l \a\2j=h-jm=l 

i-1 / .N 21 4(j-/) 

+eO) e n«-n«'. 

Z = l ^ ^ |a|2i + |i'l40-;)™=l " = ^ 

= h + l-2j 

Estimates dH]) and (ES]) state that \/h e N : 3Ch > such that HX^'^)]] < C^AT^'*. 
This also holds for h = 0, with Cq = 1. Inserted back into f l54p this shows now that 
( 153|) is true also for /c = A;^, + 1, and since fc* > 1 is arbitrary in this induction step 
while (152|) says that (153|) is true for A;^, = 1, it follows that (!53|) is true for all k E N. 
Our inductive argument also yields Rk as recursively defined for G N by 

k h~l 

Rk = E(^'^ "'"E^'^^^^^E RaiRa2^Rk-h, (57) 
/i=l 9=0 |a|2=9 

with Ro := 1 =: Co, Ci = 1/2 (09]), and C;, for /i > 2 given by coeS{X^'') at rhs(l56|). 
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To summarize, so far we have proved that the formal power series solutions for 
Dp and Hj3 satisfy 



+ \\H^''^\) < Y.RkP^''^^^^ (58) 



fc=l k=l 



with the Rk given as stated in ( 1571) and its ensuing text. We next discuss the 
convergence of rhs fj58|) . 

We begin by noting that the convergence of the series at rhsf l58l) is not altered if 
we multiply it by and add a term /^^A/" to it. Now defining ^ := /S^A/", as well as 
Ci := 1 and C2k+i '■= Rk for G N, we see that the convergence of rhs flSSl) is decided 
by the convergence of the series 

oo 

G(0:=E^2.+ie''+' (59) 

k=0 

for ,^ G C. Note that the formal power series f l59|) is just the generating function 
of the C2k+i, i-e. C2k+i = G'(2^+^)(0)/(2A; + 1)!, formally. We need to show that the 
generating function is analytic about ^ = with non-zero radius of convergence. 

With the help of the recursion relation ( [57|) we readily find that for positive ^ the 
function G{C,) is the positive inverse of the locally analytic function g ^ C, given by 

^ = ^9- '^f (60) 



which has unit derivative aX g = 0, and so its inverse ^ g is analytic in an open 
^-neighborhood of ^ = 0, with G(0) = 0. This argument already implies a finite 
radius of convergence, but without any information on its size. 

Before we actually determine the radius of convergence we note that an upper 
bound for it is readily found by discussing the function ^ defined in (1601) on 
its real interval of definition about zero. Clearly, this interval is {—g9,g,), where 



g, = y — l)/2 ^ 0.7861513775 is the positive real zero of the square root term 
in ( l60ll . Furthermore, g ^ ^ is odd, concave for positive and convex for negative 
g values, having a unique maximum at about g* = 0.4039458281 and a unique 
minimum at —g*, yielding the .^-values ^* = 0.285891853 for the maximum, and —C,* 
for the minimum; the numerical values are easily obtained using MAPLE. Clearly, 
^* is an upper bound to the radius of convergence of the power series for G{C,). 



12 



We now show that ^* is the radius of convergence. In fact, the radius of con- 
vergence is found amongst those ^ values closest to = at which the derivative 
of g vanishes (possibly asymptotically should ^ — > ^oo when g ^ oo suitably). 
Now, since any solution g = G{^) to fj60|l with G{0) = is also a solution to the alge- 
braic problem of finding the roots of the polynomial of degree 6 in g with ^-dependent 
coefficients, given by 

Pig\0 = (e - 2gYil -g'- g') - {g + g')\ (61) 

the asymptotic scenario does not occur and all one needs to do is to locate the finitely 
many zeros of the (yf-derivative of (7 ^, which reads 

^ ^ l + ?>g^ -g'-g"^ -2^1- g^-g^' ^^2) 

Clearly, a zero of r.h.s.( l62l) is also a zero of the order-12 polynomial P{g) obtained 
by multiplying the numerator of r.h.s. ( l62l) by 1 + 8(7^ — g^ — g^ + 2a/1 — g^ — g'^ , viz. 

P{9) = (1 + 3/ -g'- g'f - 4(1 - / - g'f. (63) 

By the fundamental theorem of algebra, P{g) has twelve zeros, counted in multiplic- 
ity; six of these are zeros of d^/dg given in (162 p . With the help of MAPLE one finds 
that if g^ is a zero of d^/dg given in (1^2]) . then 

^0 G ±{0.4039458281, 0.07758059914 ± il.387412147}. (64) 

Writing G^^ for r.h.s.(l60l). we now have 

|G'"^(^o)| e {0.285891853, 3.235626655}, (65) 

and the smaller of these two values is the radius of convergence of G'(^). 

Lastly, we translate the radius of convergence ^* = 0.285891853 of the series G{^) 
given in (159 p into a sufficient criterion of convergence of our formal perturbation series 
D = Dc + E^i 13^^ D^'''^ and H = Ha + E^i Z?^''^^^^- Namely, if 

(3^{\\Dc\\ + Ili^All) < 0.285891853 (66) 

then our perturbative power series solution pair for D and H converges to a classical 
solution for the stationary Maxwell-Born-Infeld field equations. 
This completes our convergence proof. 
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5 Purely electric or magnetic fields 



When J = or p = 0, then the perturbation theory for solving the stationary elec- 
tromagnetic Maxwell-Born-Infeld equations simplifies considerably. In particular, 
when applying the perturbation theory directly to the electrostatic or magnetostatic 
field equations, better convergence estimates are obtained than by simply specializing 
the electromagnetic estimates to these cases. 

5.1 The electrostatic case 

When j = and H vanishes at spatial infinity we have B = = H and the stationary 
Maxwell-Born-Infeld equations reduce to the electrostatic Maxwell-Born (-Infeld) 
equations, comprising Coulomb's law 

V ■ = 47rp , (67) 
and the stationary Faraday law 

V X E = . (68) 
The two fields E and D are linked by the aether law of Born |Borl933l IBorl969] . 

E = ^ Born's E - law. (69) 

'1 + /34|D|2 



These equations are covered in |CaKi2010] . For the convenience of the reader we 
summarize the main results in the notation of the present paper. 

Thus, inserting our perturbation theory Ansatz D = Dq+D^, with Dp = + 
^8^(2) _|_ _ ^ into (l69ll and collecting powers of /3, we find that D'^'^^ for A; G N satisfies 
the pair of linear first-order PDE 

V-D(^) = (70) 
V X = V X (71) 

where V^'^'^ is a polynomial in the D^^-* with i < k, viz. 

= -^D('=-^)^M,^ (72) 

h=l j=l \e\2j^ i=l 

h-j 
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where \^\2j and Mj have their earher assigned meaning. The pair of hnear first-order 
PDE flTOj) . flTTl) has a unique solution in Cq'" given by 

Dik) ^ pyik)^ ^ ^ 1^ (73) 

where P: Cq'" — > Cq'" projects onto the solenoidal subspace of Cg'"; see (|30|) . The 
first member of this hierarchy reads exphcitly 



Di^\x) = MDcl'Dc + V / ^ll^l}^^Si^^d\ (74) 



so /3'^( !74|) gives the leading correction to the Coulomb term D^'^^ = Dq. 

In general, the perturbation series D = ^^^q (i^^ D'^^\ with D^^^ = Dq given by 
( ITOj) and D^'^^ for e N recursively given by ( !73|) . is a formal series solution for 
the electrostatic field problem (EZ]), (EH]) with prescribed p. Following essentially the 
same strategy, in |CaKi2010] the convergence in Cq'" of the formal power series is 
proved for 

f3'^\\Dc\\ < (2^/3 - 1)^/^ ^ 0.4501964645, (75) 
in which case the perturbation series yields a classical solution D. 

5.2 The magnetostatic case 

When p = and D vanishes at spatial infinity we have D = = E and the stationary 
Maxwell-Born-Infeld equations reduce to the magnetostatic Maxwell-Born (-Infeld) 
field equations, comprising the stationary Ampere law 

V X H = Anj (76) 
and the law of the absence of sources of the magnetic induction field B, 

V ■ 5 = . (77) 
The two fields B and H are linked by the aether law of Born [ Borl933l IBorl969] , 

B = , Born's B - law. (78) 
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While these equations are not directly covered in |CaKi2010] . they can be treated 
almost verbatim to the treatment in |CaKi2010] . The minus sign under the square 
root is already covered in |CaKi2010] . and otherwise only the roles of the V- and 
Vx operators need to be exchanged. 

Thus, inserting our perturbation theory Ansatz H = + Hjj, with if^ = 
+ + . . . , into dZS]) and collecting powers of /3, we find that H^''^ for 

G N satisfies the pair of linear first-order PDE 

V X ijW = 0, (79) 
V = V-U^''\ (80) 

where U^'''' is a polynomial in the H^^'^ with i < k, viz. 

f/W = -^/7('=-^)^|M,|5^ J^ijfe-i).^fe)^ (81) 

h=l j=l le\23^ j=l 

h-3 

where again \£\2j and Mj have their earlier assigned meaning. The pair of linear 
first-order PDE (!79|) . (IHOl) has a unique solution in Cq'" given by 

i7W = Qf/W, keN (82) 

where Q : C^'" — )■ Cq'" projects onto the curl-free subspace of Cq'"; see fl3T]) . The 
first member of this hierarchy reads explicitly 



i/«(x) = i|/7Ar//A - V X / 4;FV X ,|^A^Anyj ^3^^ ^^^^ 



iFVx( 


1 

2 


I^aI 








- y 





so /3'^( l83|) gives the leading correction to the Ampere term H^^^ = Ha- 

In general, the perturbation series H = ^1^=0 with if^^^ = Ha given by 

(fTT!) and if'-'^^ for A; G N recursively given by (!82|) . is a formal series solution for the 
magnetostatic field problem (!76|) . (!77|1 with prescribed j. Following nearly verbatim 
the strategy in |CaKi2010] the convergence in Cq'" of the formal power series can be 
proved for 

I3^\\Ha\\ < (22/3-1)^/', (84) 
in which case the perturbation series yields a classical solution H . 
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6 Summary and outlook 



We have constructed a convergent perturbative power series which solves the sta- 
tionary Maxwell-Born-Infeld equations with regular sources p and j in Cg'" as long 
as l3'^{\\Dc\\ + \\Ha\\) < r = 0.28 (approximately); here, || . || is the Cg'" norm, Dq 
the Coulomb field and the Ampere field. In the absence of either Dc or Hp^ the 
corresponding perturbative power series for the remaining field even converges when 
/3W < ^ 0.45, where M is the Cq'" norm of i^A, resp. Dc] cf. |CaKi2010] . 

While our solution method does not cover the conceptually important case of delta 
function sources for which Born conceived of this nonlinear field theory |Borl933| 
IBorl937] IBoInl933] , it nevertheless allows one to gain some relevant insight into a 
burning open question: reliable hounds on the size of Bom's aether constant (3. 

In a purely electrostatic calculation Born and Infeld [BoInl933] computed /3 = 
/3b ~ 1.236as (in the units of |Kie2004al [Kie2004bj . where ag ^ 1/137.036 is Som- 
merfeld's fine-structure constant; so, (3b ~ 0.00902. Subsequently, trying to take the 
magnetic moment of the electron into account. Born and Schrodinger |BoSchrl935] 
came up with a purely magnetostatic (very rough) estimate of (3 = /3bs ~ 4.74/3b. 
Subsequently Rao |Raol936] undertook a more detailed electromagnetic study, fol- 
lowing Born's suggestion to treat the electron not as a point but as a charged, rotating 
(one-dimensional) ring. However, all these plausible estimates are not based on an 
experimentally accessible process, involving as they do solely a single "stationary 
electron solution" in an otherwise empty space. Various attempts have been made 
since to compute effects on atomic spectra induced by the Born-Infeld nonlinear- 
ity to estimate the range of viable (3 values (see |CaKi2006] and references therein; 
see also |FrGa20lO] for a recent study), but all these methods are based on various 
as-of-yet uncontrolled approximations and have produced conflicting results. 

Our solution method now paves the way for a systematic study based on controlled 
approximations. For example, if we replace the Dirac delta function of each point 
charge by a C°'^ function supported on a ball of radius equal to the Compton wave- 
length of the electron (which sets the unit of length in |Kie2004at IKie2004b] ). then 
\\Dc\\ ~ 3q;s, so our electrostatic series will converge for (3 ^ 0.3874/ ^/as ~ 4.5348, 
which range of /3 values safely includes the value proposed by Born. In particular, 
our perturbation series should converge sufficiently rapidly so that only a few terms 
need to be computed; see the appendix. Of course, to not compare apples with or- 
anges, the resulting spectral data computed for the Schrodinger operators obtained 
with the help of these systematically computed D = Dc + Dp flelds for smeared 
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out charges should be compared with those computed with the Schrodinger operator 
with Coulomb potential for the same "smeared-out" charges, not point charges. 

Similarly, we may also try to estimate the viable range of /3 values for an electron 
model with magnetic moment, although comparison with the purely magnetostatic 
calculation by Born-Schrodinger will at best be an academic excercise because such 
a model (presumably) can't produce the Hydrogen spectrum in the /3 — )■ limit. 
Instead one should base the assessment on the computation of a stationary electro- 
magnetic solution with both p and j sources supported on a "smeared-out" ring, in 
the spirit of Rao's study |Raol936] . However, one possible caveat is that the {3 value 
may get uncomfortably close to the radius of convergence of the electromagnetic 
series, so that more terms may need to be calculated than is practically feasible. 

We close the main body of this paper with an outlook on desired generalizations 
of the perturbative solution strategy pursued here for the stationary Maxwell-Born- 
Infeld field equations in three-dimensional Euclidean space. The generalization to the 
stationary model in higher- dimensions, which in recent years has become of interest 
to the high energy theory community (see, e.g., |Gibl998t ILiYa2003] ) . should be quite 
straightforward (cp. the electrostatic version in [CaKi2010] ). Less straightforward, 
but very desirable, is the generalization to the genuinely time-dependent setting; 
see |Serl988l ISer2004l IBren2004l |Spe2010| for recent rigorous results, |Boil970j for 
an early pioneering study, and |BiBil983] for a delightful general discussion. The 
main problem to be overcome to make the perturbative strategy work in the dy- 
namical setting is the identification of a suitable Banach algebra. Lastly, and most 
importantly from the Born-Infeld perspective [Borl933l IBoInl933] . it is desirable 
to have the method which can handle delta function sources both in fiat space (see 
[Borl933l IBoInl933l IEckl986j and |Hopl994i IGibl998j for some explicit solutions 
with a single, respectively infinitely many point charges) and in curved spacetime 
(see |Hofl935al IHofl935bl IGSP19841 IDeml986l IBret2001l ITaZa2010] for explicit so- 
lutions with a single point charge). The technique presented in this paper can serve 
to produce seed field-configurations for growing a solution for delta function sources 
in a suitable limiting process. Progress in this direction will be reported elsewhere. 
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opinions expressed in this paper are those of the author and do not necessarily reflect 
those of the NSF. I thank Shadi Tahvildar-Zadeh and Holly Carley for their com- 
ments. Special thanks are due to Doron Zeilberger and Shalosh B. Ekhad |Ekh2050] 
for their wizzardry help with the recursion relation for the coefficients Rk- 
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Appendix^ 



Practically only a small number of terms of the perturbation series for D and H 
can be computed. Here we supply some details about the recursion coefficients Rk 
which are needed in the error estimates. Given the input value i?o = 1, the author 
computed Ri,R2,R3 by hand from the recursion formula given in the main text, 
then used MAPLE to compute 14 terms using the iteration of the fixed point map 
given by the implicit definition of G{^) in a direct, unpolished manner (after which 
MAPLE reached its capacity). The more sophisticatedly iterating Shalosh B. Ekhad 
knows 303 terms of the sequence {-RfejfcgNujo}; here are the first 30 of them: 

i?o = 1 
3 



Ri 

R2 

R3 
R4 

R^ 
R& 
R7 
Rg 

Rq 

Rio 
Rii 
R12 



2 

65 

~8 
943 

62689 

128 
1128197 

256 
42790845 

1024 
842157399 

2048 
136312116961 

32768 
2817640708457 

65536 
118490151386655 

262144 
2526390089218393 

524288 
217977129447815405 

4194304 



^5) 
^6) 

^7) 



89) 
90) 
91) 
92) 
93) 
94) 
95) 
96) 
97) 



^Disclaimer: Everything in this appendix could be made rigorous, but this here is only semi- 
rigorous. The whole process took Ekhad 200.542 seconds. 



19 



_ 4748017126294329161 

^ 8388608 ^^^^ 

208581111836961908919 

- 33554432 ^ ' 

4614991020517094410279 , 

Ri5 = (100 

^ 67108864 ^ ^ 

1644116252728526666074977 

- 2147483648 

_ 36812969231234813601419473 

^"^ ~ 4294967296 

1656740336870818323274233515 

- 17179869184 

37445969415289365495538129125 , 

Rig = 104 

^ 34359738368 ^ ' 

3398982473269915594232889691503 

- 274877906944 

77410530113072758320538102052283 
R21 = (106) 

549755813888 ^ ^ 

_ 3537571318060518004220386126288923 
2199023255552 

81073698522885193141789999978889377 , , 

R23 = 108 

4398046511104 ^ ^ 

14905122955618940253574385037312323437 

R24 = (109) 

70368744177664 ^ ^ 

343396823100629008332240991968973221221 , , 

R25 = 110 

140737488355328 ^ ^ 

15859792436328056179243840618808531803115 , , 

R26 = 111 

^ 562949953421312 ^ ^ 

367032135637139188746851720898633881475805 

~ 1125899906842624 ^ ^ 

_ 34043872901750574940303463222092806831775365 

~ 9007199254740992 

790900804799227563339255535454710815651794705 

R29 = ■ (114) 

^ 18014398509481984 ^ ^ 

For many practical purposes the above list should be fully sufficient. However, 
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should more terms be needed it is more efficient to switch to an asymptotic repre- 
sentation. Unfortunately, the recursion formula found in this paper is not practically 
useful for this purpose, for it quickly grows out of control because it involves all previ- 
ous Re values {i < k) to compute the k-th. one. And the fix point iteration produces 
the Rk without yielding control of their asymptotics. Fortunately, we know from 
general theorems that the associated sequence of integers {Sk+i '■= '^'^'^ Rk} keNu{o} 
satisfies a linear recurrence of fixed finite order with coefficients which are fixed 
polynomials in n = A; + 1 G N. The beginning of this integer sequence reads 
{1, 6, 130, 3772, 125378, ...}. After a few minutes of interaction with Doron Zeilberger, 
Shalosh B. Ekhad found: given {S*!, S'g}, the ensuing Sn values are obtained as 

Sn+1 = / Ph{n)Sn^h, (115) 

where each Ph is a polynomial of degree 9 in n, with integer coefficients. The coeffi- 
cients of these polynomials are insanely large, with more than 100 digits, so that it 
would be pointless to list them here; they are available from Dr.Z.'s computer. 

The asymptotic behavior of {Sfc+ijfegNuio} can now be extracted from the linear 
recurrence. The leading term in the asymptotics is determined by the degree-8 
polynomial in the complex variable z formed by collecting the coefficients of highest 
power (i.e. n^) of the polynomials Ph{n) in (IllSp . each multiplied by the associated 
h-th power of z, thus (symbolically) 

where P^\n) is the 9-th derivative of Ph{n), with n formally treated as real variable. 
As 77, — )■ cxD, the ratio Sn+i/Sn converges to the largest root of flll6p . which Ekhad 
easily finds as ~ 48.9391511596 (truncated after 10 digits). Dividing by 4 yields 

lim Rk+i/Rk = g ~ 12.23478778990. (117) 

fc— >oo 

As many digits of the ratio q are practically reached only for /c ^ 30; however, 
empirically, R^+i/Rk is monotonic increasing, and so 

Rk < q^. (118) 

This suffices to obtain the relevant error estimates for the truncation of the series. 

Incidentally, by the ratio criterion, flll7p yields the radius of convergence of G(0 
as ^* = l/y^ ~ 0.285891853, reproducing what we found earlier quasi-algebraically. 
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